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ABSTRACT 

Stars form in dense cores of molecular clouds that are observed to be signif- 
icantly magnetized. In the simplest case of a laminar (non-turbulent) core with 
the magnetic field aligned with the rotation axis, both analytic considerations 
and numerical simulations have shown that the formation of a large, 10^ AU- 
scale, rotationally supported protostellar disk is suppressed by magnetic braking 
in the ideal MHD limit for a realistic level of core magnetization. This theoreti- 
cal difficulty in forming protostellar disks is termed ^^magn etic braking catastro- 



phe. " A possib le resolution to this problem, proposed by iHennebelle fc Ciardi 



and lJoos et al.l . is that misalignment between the magnetic field and rotation axis 
may weaken the magnetic braking enough to enable disk formation. We evaluate 
this possibilit y quantitat ively through numerical simulations. We confirm the 



basic result of iJoos et al.l that the misalignment is indeed conducive to disk for- 
mation. In relatively weakly magnetized cores with dimensionless mass-to-fiux 
ratio > 5, it enabled the formation of rotationally supported disks that would 
otherwise be suppressed if the magnetic field and rotation axis are aligned. For 
more strongly magnetized cores, disk formation remains suppressed, however, 
even for the maximum tilt angle of 90°. If dense cores are as strongly magnetized 
as indicated by OH Zeeman observations (with a mean dimensionless mass-to- 
fiux ratio ^ 2), it would be difiicult for the misalignment alone to enable disk 
formation in the majority of them. We conclude that, while beneficial to disk 
formation, especially for the relatively weak field case, the misalignment does not 
completely solve the problem of catastrophic magnetic braking in general. 
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Introduction 



Star and planet formation are connected through disks. Disk formation, long thought 
to be a trivial c onsequence of angu lar momentum conservation during core collapse and star 
formation (e.g. . lBodenheimerlll995l ) . turned out to be much more complicated than originally 
envisioned. The complication comes fro m magnetic field s, which are observed in dense, star- 
forming, cores of molecular clouds (see ICrutcherl l2012l for a recent review). The field can 
strongly affect the angular momentum evolution of core collapse and disk formation through 
magnetic braking. 

There have been a number of studies aiming at quantifying the effects of magnetic 
field on disk formation. In the ideal MHD limit, both analytic considerations and nu- 
merical simulations have shown that the formation of a rotationally supported disk (RSD 
hereafter) is suppressed by a realistic magnetic field (corresponding to a dimensionless mass- 
to-fiux ratio of A ^ a few; iTroland fc Crutcherll2008l ) during the protostellar mass accretion 



rotation axis ( 


Allen et al. 


2oo: 




Galli et al. 




2006 




Price & Bate 


2007 




Hennebelle & Fromang 


20081: ] 


3aDD & Basu 


2010 


Seifried et al.l 


2011 




2012). The suppression of RSDs bv excessive 


magnetic braking is termed 



Mellon & Li 




2008; 


Santos- Lima et al. 



catastrophe" in star formation. 

Rotationally supporte d disks are routinely obs erved, however, around evolved Class II 



young stellar objects (see 



W 



Class I /e.g.. |j0rgensen et al 



Lee 



lliaras fc Cieza 
2OO9I : 



2011 



2011 



for a review), an d increasingly around 



Takakuwa et al.ll2012l ) and even one Class 



source JXcbin et aljl2012h . When and how such disks form in view of the magnetic braking 



catastrophe is unclear. The current attempts to overcome the catastrophic braking fall into 
three categories: (1) non-ideal MHD effects, including ambipolar diffusion, Ohmic dissipa- 
tion and Hall effect, (2) misalignment between magnetic and rotation axes, and (3) turbu- 
lence. Ambipolar diffusion does not appear to weaken the braking enough to enable large- 
scale RSD formation under r e alistic conditio ns ( IKrasnopolsky fc Konigll l2002l : Mellon fc Li 
2OO9I : iDuffin fc Pudritzl l2009l : iLi et al.l I2OIII ) . Ohmic dissipation can prod u ce small, AU- 
scale , RSD in the early protostellar accreti on phase ( IMachida et al.l I2OIOI : iDapp fc Basu 



2OIOI : iDapp et al 



2012: 



Tomida et al. II2013I ). Larger, 10^ AU-scale RSDs can be produced 
if the resistivity or t he Hall coefficient of the dense core is n iuch larger than the classi cal 
(microscop i c) val ue ( IKrasnopolsky et al.l 1201 Ol . 1201 ll : see also braiding fc Wardlel l2012al lbl) . 
Joos et al.l ( 120 12I ) explored th e effects of tilting the magnetic field away from the rotation 
axis on disk formation (see also lMachida et al.lbood : IPrice fc BatlboO?! : iHennebelle fc Ciardil 
2OO9I ) . They concluded that Keplerian disks can form for a mass-to-fiux ratio A as low as 
3, as long as the tilt angle is close to 90° (see their Fig. 14). The effects of turbulence were 
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explored by ISantos-Lima et al.l ( 120121 l2013l ) , who concluded that a strong enough turbu- 
lence can induce eno ugh magnetic diffusion t o enable the formation of a 10^ AU-scale RSD. 
Seifried et al.l ( 120121 ) and iMyers et al.l ( 120121 ) considered supersonically turbulent massive 
cores. They found rotationally dominated disks around low-mass stars, although in both 
cases the turbulence-induced rotation is misaligned with the initial magnetic field b y a large 
angle, which may have contributed to the disk formation (see also lJoos et al.ll2013l ). 



The goal of this paper is to revisit the role of magnetic field-rotation misalignment in disk 
formation. The misalignment is expected if the angular momenta of dense cores are generated 
through turbulent motions (e.g.. iBurkert fc Bodenheimerll2000l : iMyers et al.ll2012l ). It is also 
inferred from the misalignment between the field direction traced by polarized du st emission 



and the outfiow axis, which is taken a s a proxy f or the direction of rotation ( IHull et al. 



20121). Indeed, in the CARMA sample of IHuU et al.l . the distribution of the angle between 
the magnetic field and jet/rotation axis is consistent with being random. If true, it would 
indicate that in half of the sources the two axes are misaligned by a large angle of 9q > 60° 
(see however IChapman et al.ll2013l and discussion in §5.21) . Such a large misalignment would 
be enoug h to allow disk forraation in dense cor es magnetized to a realistic level (with A 
of a few; iTroland fc Crutcherl l2008l) according to Ijoos et al.l (120121 ). If the alignment angle 
^0 is indeed random and IJoos et al.r s conclusions are generally true, the magnetic braking 
catas trophe woul d be largely solved. Given their far-reaching implications, it is prudent to 
check I Joos et al.r s conclusions, using a different numerical code. It is the task of this paper. 



We carry out numerical experiments of disk formation in dense cores with misaligned 
magnetic and rotation axes using non-ideal MHD code Zeus-TW that includes self-gravity. 
We find that a large misalignment angle does indeed enable the formation of RSDs in weakly 
magnetized dense cores with dimensionless mass-to-fiux ratios > 5, but not in dense cores 
magnetized to higher, more typical levels. Our conclusion is that while the misalignment 
helps with disk formation, especially in relatively weakly magnetized cores, it may not provide 
a complete resolution to the magnetic braking catastrophe by itself. 

The rest of the paper is organized as follows. In ^ we describe the model setup. The 
numerical results are described in §3] and §11 We compare our results to those of I Joos et al. 
and discuss their implications in §5] and conclude with a short summary in §5.31 



2. Problem Setup 



We follow iLi et al.l (120111 ) and iKrasnopolskv et al.l (120121 ) and start our simulations 
from a uniform, spherical core of 1 Mq and radius 10^'' cm in a spherical coordinate system 
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(r, 0,0). The initial density po = 4.77 x 10~^^gcm~^ corresponds to a molecular hydrogen 
number density of lO^cm"^. We adopt an isothermal equation of state with a sound speed 
a = 0.2kms~^ below a critical density Pc = 10~^^gcm~^, and a polytropic equation of state 
p oc p^/^ above it. At the beginning of the simulation, we impose a solid-body rotation 
of angular speed Qq = 10~^^s~^ on the core, with axis along the north pole (9 = 0). It 
corresponds to a ratio of rotational to gravitational binding e nergy of 0.025, which is typical 



of the values inferred for NH3 cores ( IGoodman et al.lll993l ). The initial magnetic field is 



uniform, tilting away from the rotation axis by an angle ^o- We consider three values for the 
initial field: Bq = 10.6, 21.3 and 35.4 /xG, corresponding to dimensionless mass-to-fiux ratio, 
in units of {27tG^/^)-\ A = 9.72, 4.86 and 2.92, respectively, for the core as a whole. The 
mass-to-fiux ratio for the central fiux tube Ac is higher than the global value A by 50%, so that 
Ac = 14.6, 7.29 and 4.38 for the three cases respectively. The effective mass-to-fiux ratio Agff 
should lie betwee n these two limits. If the star formation eflSciency per core is ^ 1/3 (e.g.. 



Alves et al.l 120071 ). then one way to estimate Agff is to consider the (cylindrical) magnetic 
fiux surface that encloses 1/3 of the core mass, which yields Aeff = 1.41 A, corresponding to 
13.7, 6.85, and 4.12 for the three cases respectively; the fraction 1/3 is also not far from 
the typical fraction of core mass that has accreted onto the central object at the end of our 
simulations (see Table 1). For the tilt angle, we also consider three values: 9q = 0°, 45° and 
90°. The ^0 = 0° corresponds to the aligned case, with the magnetic field and rotation axis 
both along the z-axis (^ = 0). The 9q = 90° corresponds to the orthogonal case, with the 
magnetic field along the x-axis {9 = 90°, = 0). Models with these nine combinations of 
parameters are listed in Table 1; additional models are discussed below. 



As in iKrasnopolsky et al.l (120121 ). we choose a non-uniform grid of 96 x 64 x 60. In 
the radial direction, the inner and outer boundaries are located at r = 10^^ and 10^ ''cm, 
respectively. The radial cell size is smallest near the inner boundary (5 x 10^^ cm or ^ 
0.3 AU). It increases outward by a constant factor ^ 1.08 between adjacent cells. In the 
polar direction, we choose a relatively large cell size (7.5°) near the polar axes, to prevent the 
azimuthal cell size from becoming prohibitively small; it decreases smoothly to a minimum of 
^ 0.63° near the equator, where rotationally supported disks may form. The grid is uniform 
in the azimuthal direction. 

The boundary conditions in the azimuthal direction are periodic. In the radial direction, 
we impose the standard outfiow boundary conditions. Material leaving the inner radial 
boundary is collected as a point mass (protostar) at the center. It acts on the matter in the 
computational domain through gravity. On the polar axes, the boundary condition is chosen 
to be refiective. Although this is not strictly valid, we expect its effect to be limited to a 
small region near the axis. 
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Table 1. Models 



Model 


A'* 


eiT 




1] (lO^'^cm^s-^) 


(Mo) 


RSD'^ 


A 


9.72 


13.7 


0° 


1 


0.24 


No 


B 


4.86 


6.85 


0° 


1 


0.22 


No 


C 


2.92 


4.12 


0° 


1 


0.33 


No 


D 


9.72 


13.7 


45° 


1 


0.21 


Yes/ Porous 


E 


4.86 


6.85 


45° 


1 


0.35 


No 


F 


2.92 


4.12 


45° 


1 


0.27 


No 


G 


9.72 


13.7 


90° 


1 


0.38 


Yes/ Robust 


H 


4.86 


6.85 


90° 


1 


0.46 


Yes/ Porous 


I 


2.92 


4.12 


90° 


1 


0.47 


No 


M 


9.72 


13.7 


90° 





0.10 


Yes/ Robust 


N 


9.72 


13.7 


90° 


0.1 


0.26 


Yes/ Robust 


P 


4.03 


5.68 


90° 


1 


0.25 


Yes/ Porous 


Q 


3.44 


4.85 


90° 


1 


0.14 


No 



Note. — a). The average dimensionless mass-to-flux ratio for the core as a whole; b). The 
eflFective mass-to-flux ratio for the central 1/3 of the core mass (see ^ for discussion); c). 
Mass of the central object when the simulation is stopped; d). "Robust" disks are persistent, 
rotationally supported structures that rarely display large deviations from smooth Keplerian 
motions, whereas "porous" disks are highly active, rotationally dominated structures with 
large distortions and may occasionally be completely disrupted. 
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We initially intended to carry out simulations in the id eal MHD l i mit, s o that they 
can be compared more directly with other work, especially Ijoos et al.l (120121 ). However, 
ideal MHD simulations tend to produce numerical "hot zones" that force the calculation to 
stop early in the proto stellar mass accretion phase, a tendency we no ted in our previous 
2D dMellon fc Lil boosi ) and 3D simulations krasnopolskv et al.l [2OI2I ). To lengthen the 
simulation, we include a small, spatially uniform resistivity r] = lO-^^cm^s"^. We have 
verified that, in the particular case Model G (A = 9.72 and 6^0 = 90°), this resistivity changes 
the flow structure little compared to either the ideal MHD Model M (before the latter stops) 
or Model N, where the resistivity is reduced by a factor 10, to 10^ 



3. Weak-Field Case: Disk Formation Enabled by Field-Rotation Misalignment 

3.1. Equatorial Pseudodisk vs Magnetically Induced Curtain 

To illustrate the effect of the misalignment between the magnetic field direction and 
rotation axis, we first consider an extreme case where the magnetic field is rather weak 
(with a mass-to-fiux ratio A = 9.72 for the core as a whole and Aeff = 13.7 for the inner 
1/3 of the core mass). In this case, a well-formed rotationally supported disk is present 
in the orthogonal case with = 90° (Model G in Table 1). Such a disk is absent in the 
aligned case (with ^0 = 0°, Model A). The contrast is illustrated in Fig. [U where we plot 
snapshots of the aligned and orthogonal cases at a representative time t = 3.9 x 10^^ s, when 
a central mass of 0.11 and 0.12 M©, respectively, has formed. The fiow structures in the 
two cases are very different in both the equatorial (panels [a] and [b]) and meridian (panels 
[c] and [d]) plane. In the equatorial plane, the aligned case has a relatively large (with 
radius ^ 10^^ cm) over-dense region where material spirals rapidly inward. On the (smaller) 
scale of 10^^ cm, the structure is dominated by expanding, low-density lobes; they are the 
decoupling enable d m agnetic structures (DBMS f or short) that have been studied in detail by 



Zhao et al.l ( 120111 ) and lKrasnopolsky et al.l ( 120121 ). No rotationally supported disk is evident. 



The equatorial structure on the 10^^ cm scale in the orthogonal case is dominated by a pair 
of spirals instead. The spirals merge, on the 10^^ cm scale, into a more or less continuous, 
rapidly rotating structure — a rotationally supported disk. Clearly, the accretion fiow in the 
orthogonal case was able to retain more angular momentum than in the aligned case. Why 
is this the case? 

A clue comes from the meridian view of the two cases (panels [c] and [d] of Fig. [1]). 
In the aligned case, there is a strong bipolar outfiow extending beyond 3 x 10^^ cm at the 
relatively early time shown. The outfiow forces most of the infalli nR material to acc rete 



through a fiattened equatorial structure — an over-dense pseudodisk (iGalli fc Shu 



1993; see 



-MO^^ -5«10^^ 5«10^^ MO^^ -MO^^ -5«10^^ 5«10^^ MO' 




-3«10^^ -2«10^^ -MO^^ MO^^ 2«10^^ 3«10^^ -3«10^^ -2«10^^ -MO^^ MO^^ 2«10^^ 3«10' 



Fig. 1. — Snapshots of the logarithm of density (color map, in units of gcm~^) and velocity 
field (white arrows) for the weakly magnetized core of A = 9.72 at a representative time 
t = 3.9 X 10^^ s. The left (right) panels are for the aligned (orthogonal) case, and the top 
(bottom) panels are for the equatorial x-y (meridian y-z) plane. Note the powerful bipolar 
outflow in panel (c) driven by the pseudodisk in panel (a). The lack of powerful outflow in 
panel (d) is indicative of a weaker magnetic braking, which is consistent with the presence 
of a rotationally supported disk in panel (b). The length unit is cm. The scale in (c) and 
(d) is three times larger than that in (a) and (b). 
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panel [a] for a face-on view of the pseudodisk, noting the difference in scale between panel 
[c] and [a]). It is the winding of the magnetic field lines by the rotating material in the 
pseudodisk that drives the bipolar outfiow in the first place. The wound-up field lines act 
back on the pseudodisk material, braking its rotation. It is the efficient magnetic braking in 
the pseudodisk that makes it diflScult for rotationally supported disks to form in the aligned 
case. 

The prominent bipolar outffow indic ative of eflScient mag n etic b raking is absent in the 



orthogonal case, as was emphasized by ICiardi fc Hennebelld ( 120101 ). It is replaced by a 



much smaller, shell-like structure inside which the 10^^ cm-scale rotationally supported disk 
is encased (panel [d]). To understand this difference in ffow structure pictorially, we plot in 
Fig. [2] the three-dimensional structure of the magnetic field lines on the scale of 1000 AU (or 
1.5 X 10-^^ cm), which is 50% larger than the size of panels (a) and (b) of Fig.dJ but half of that 
of panels (c) and (d). Clearly, in the aligned case, the relatively weak initial magnetic field 
(corresponding to A ?^ 10) has been wound up many turns by the material in the equatorial 
pseudodisk, building up a magnetic pressure in the equatorial region that is released along 
the polar directions (see the first panel of Fig. [2]). The magnetic pressure gradient drives a 



Tomisaka ( 


1998) and 


Allen et al. 


( 


2003) 



In contrast, in the orthogonal case, the equatorial region is no longer the region of the 
magnetically induced pseudodisk. In the absence of rotation (along the z-axis), the dense 
core material would preferentially contract along the field lines (that are initially along the 
X-axis) to form a dense sheet in the y-z plane that passes through the origin. The twisting of 
this sheet by rotation along the z-axis produces two curved "curtains" that spiral smoothly 
into the disk at small distances, somewhat analogous in shape to two snail-shells (see the 
second panel of Fig. [J]). 

The snail-shaped dense curtain in the orthogonal case naturally explains the morphology 
of the density maps shown in panels (b) and (d) of Fig. [H First, the two prominent spiral arms 
in the panel (a) are simply the equatorial {x-y) cut of the curved curtains. An interesting 
feature of the spirals (and the snail-shaped dense curtain as a whole) is that they are the 
region where the magnetic field lines change directions sharply. This is illustrated in Fig. 
[3l which is similar to panel (b) of Fig. [U except that the magnetic vectors (rather than 
velocity vectors) are plotted on top of the density map. Clearly, the spirals separate the 
field lines rotating counter clock-wise (lower-right part of the figure) from those rotating 
close- wise (upper-left). The sharp kink is analogous to the well-known field line kink across 
the equatorial pseudodisk in the aligned case, where the radial component of the magnetic 
field changes direction. It supports our interpretation of the spirals and, by extension, the 
curtain as a magnetically induced feature, as is the case of pseudodisk. In other words, the 




Fig. 2. — 3D view of representative magnetic field lines and isodensity surfaces at p = 10~^^ 
(red) and 10~^^gcm~^ (blue) for the aligned Model A and orthogonal Model G. For clarity, 
only field lines originated from the bottom x-y and left y-z plane are plotted, respectively. 
Note that the magnetically induced equatorial pseudodisk in the aligned case is warped by 
rotation into two snail-shaped curtains that spiral inward to form a disk (blue surface) in 
the orthogonal case. The length is in units of AU. 
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spirals are not produced by gravitational instability in a rotation ally supported stru cture; 



they are "pseudospirals" in the same sense as the "pseudodisks" of iGalli fc Shul ( 119931 ) . The 
field line kinks are also evident across the dense curtain in the 3D structure shown in the 
second panel of Fig. [21 

The 3D topology of the magnetic field and the dense structures that it induces lie at the 
heart of the difference in the magnetic braking efficiency between the aligned and orthogonal 
case. In particular, a fiattened, rotating, equatorial pseudodisk threaded by an ordered 
magnetic field with an appreciable vertical component (along the rotation axis) is more 
conducive to driving an outfiow than a warped curtain with a magnetic field predominantly 
tangential to its surface. The outfiow plays a key role in angular momentum removal and 
the suppression of rotationally supported disks, as we demonstrate next. 



3.2. Torque Analysis 



To quantify the outfiow effect, we follow IZhao fc Lil ( 120131 : see also lJoos et al.ll2012l ) and 
compare the rates of angular momentum change inside a finite volume V through its surface 
S due to infall and outfiow to that due to magnetic torque. The total magnetic torque 
relative to the origin (from which a radius vector r is defined) is 



Nm = ^ j[rx{{VxB)xB)]dV, 



where the integration is over the volume V. Typically, the magnetic torque comes mainly 
from the magnetic tension rather than pressure force. The dominant magnetic tension term 
can be simplified to a surface integral ( IMatsumoto fc Tomisakall2004l ) 



N* = ^ y"(rxB)(B.dS), 



(2) 



over the surface S of the volume. This volume-integrated magnetic torque is to be compared 
with the rate of angular momentum advected into the volume through fluid motion, 



J p(rxv)(vdS), 



(3) 



which will be referred to as the advective torque below. 



Since the initial angular momentum of the dense core is along the z-axis, we will be 
mainly concerned with the z-component of the magnetic and advective torque which, for a 
spherical volume inside radius r, are given by 

1 

4:71 



(4) 



Fig. 3. — Same as panel (b) of Fig. [1] but with magnetic unit vectors instead of velocity 
vectors plotted. Note that the field lines kink sharply in the spirals. The kinks demonstrate 
that the spirals (and the curtain that contains them) are the counterparts to the equatorial 
pseudodisk in the aligned case. 



and 



I 



pZUV(f)VrdS. 



(5) 



The advective torque consists of two parts: the rates of angular momentum advected into 
and out of the sphere by infaU and outflow respectively: 



An example of the magnetic and advective torques is shown in Fig. [H The torques are 
evaluated on spherical surfaces of different radii, at the representative time t = 3.9 x 10^^ s. 
For the aligned case, the net torque close to the central object is nearly zero up to a radius of 
^ 2 X 10^^ cm, indicating that the angular momentum advected inward is nearly completely 
removed by magnetic braking there. At larger distances, between ~ 2 x 10^^ and ^ 10^^ cm, 
the net torque Nt^z + A^a,z is negative, indicating that the angular momentum of the material 
inside a sphere of radius in this range decreases with time. This is in sharp contrast with 
the orthogonal case, where the net torque is positive in that radial range, with the angular 
momentum there increasing (rather than decreasing) with time. One may think that the 
difference is mainly due to a significantly larger magnetic torque Nt^z in the aligned case 
than in the orthogonal case. Although this is typically the case at early times, the magnetic 
torques in the two cases become comparable at later times (see the lowest solid lines in 
the two panels of Fig. [H a movie of the torques is available on request from the authors). 
A bigger difference comes from the total (or net) angular momentum Na^z = A^a^^ + A^^^* 
advected inward, which is substantially smaller in the aligned case than in the orthogonal 
case (see the uppermost solid lines in the two panels). The main reason for the difference 
is that a good fraction of the angular momentum advected inward by infall A^^^ is advected 
back out by outflow A^^* in the former, but not the latter. This is helped by the fact that 
A^^^ is somewhat smaller in the aligned case to begin with (compare the dotted lines in the 
two panels). The lack of appreciable outward advection of angular momentum by outflow, 
which is itself a product of field- winding and magnetic braking in the aligned case, appears to 
be the main reason for the orthogonal case to retain more angular momentum at small radii 
and form a rotationally supported disk in this particular case of relatively weak magnetic 
field. 




(6) 



and 




(7) 



The formation of a rotationally supported disk can be seen most clearly in Fig. [5l where 
we plot the infall and rotation speed, as well as specific angular momentum as a function of 
radius along 4 (±x and ±y) directions in the equatorial plane. In the orthogonal case, the 
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Fig. 4. — Magnetic (black solid line, A^t,^) and total advective (magenta, Na^z) torque acting 
on spheres of different radii for the aligned (left panel) and orthogonal (right) case at the 
representative time t = 3.9 x 10^^ s. Also plotted are the net torque (red, Nt^z + ^a.z)^ 
and the contributions to the advective torque by infall (dotted, N^J^^) and outflow (dashed, 
N^^^). Note that the net torque is negative between ^ 2 x 10^^ and ^ 10^^ cm for the aligned 
case but positive for the orthogonal case, mainly because the outflow removes more angular 
momentum in the former than in the latter. 
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infall and rotation speeds display the two tell-tale signs of rotationally supported disks: (1) 
a slow, subsonic (although nonzero) infall speed much smaller than the free fall value, and 
(2) a much faster rotation speed close to the Keplerian value inside a radius of ^ 100 AU. 
The absence of a rotationally supported disk in the aligned case is just as obvious. It has 
a rotation speed well below the Keplerian value and an infall speed close to the free fall 
value, especially at small radii up to ^ 100 AU. This corresponds to the region dominated 
by the low-density, strongly magnetized, expanding lobes (i.e., DEMS; see panel [a] of Fig[l]) 
where the angular momentum is almost completely removed by a combination of magnetic 
torque and outflow (see the third panel in Fig. E]). Also evident from the panel is that the 
speciflc angular momentum of the equatorial inflow drops signiflcantly twice: near ^ 10^^ 
and ^ 10^^ cm respectively. The former corresponds to the DEMS-dominated region, and 
the latter the pseudodisk (see panel [a] of Fig. [1]). The relatively slow infall inside the 
pseudodisk allows more time for magnetic braking to remove angular momentum. It is the 
pseudodisk (and its associated outflow) working in tandem with the DEMS that suppresses 
the formation of a rotationally supported disk in the aligned case. Interestingly, there is a 
bump near ^ 10^^ cm for the speciflc angular momentum of the orthogonal case, indicating 
that the angular momentum in the equatorial plane is transported radially outward along 
the spiraling fleld lines from small to large distances (see Fig. [3]). 

The spiraling equatorial fleld lines in the orthogonal Model G have an interesting prop- 
erty: they consist of two strains of opposite polarity. As the strains get wound up more 
and more tightly by rotation at smaller and smaller radii, fleld lines of opposite polarity 
are pressed closer and closer together, creating a situation that is conducive to reconnec- 
tion, either of physical or numerical origin (see the flrst panel of Fig. [6]). Model G contains 
a small but flnite resistivity {rj = lO^'' cm^ s~^). It does not appear to be responsible for 
the formation and survival of the Keplerian disk, because a similar disk is also formed at 
the same (relatively early) time for a smaller resistivity of rj = 10^^ cm^ s~^ (Model N) and 
even without any explicit resistivity (Model M). Numerical resistivity may have played a 
role here, but it is diflicult to quantify at the moment. In any case, the magnetic fleld on 
the Keplerian disk appears to be rather weak, as can be seen from the second panel of Fig. 
[6l where the plasma /3 is plotted along 4 (±x and ±y) directions in the equatorial plane. 
On the Keplerian disk in the orthogonal Model G (inside ~ 100 AU), /3 is of order 10^ or 
more, indicating that there is more matter accumulating in the disk than magnetic fleld, 
either because the matter slides along the fleld lines into the disk (increasing density but 
not the fleld strength) or because of numerical reconnection that weakens the fleld, or both. 
This situation is drastically different from the aligned case, where the inner 100 AU region 
is heavily influenced by the magnetically-dominated low-density lobes. 
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4. Moderately Strong Field Cases: Difficulty with Disk Formation 

We have seen from the preceding section that, in the weakly magnetized case of A = 9.72, 
the 10^ AU-scale inner part of the protosteUar accretion flow is dominated by two very 
different types of structures: a weakly magnetized, dense, rotationally supported disk (RSD) 
in the orthogonal case (^o = 90°, Model G) and magnetically dominated, low-density lobes or 
DEMS in the aligned case (^o = 0°, Model A), at least at the relatively early time discussed 
in §31 when the central mass reaches ^ 0.1 Mq. This dichotomy persists to later times for 
these two models and for other models as well, as illustrated by Fig. [71 where we plot Models 
A-I at a time when the central mass reaches 0.2 Mq. It is clear that the RSD for Model G 
becomes even more prominent at the later time, although a small magnetically dominated, 
low-density lobe is evident close to the center of the disk: it is a trapped DEMS that is too 
weak to disrupt the disk. In this case, the identiflcation of a robust RSD is secure, at even 
later times (up to the end of the simulation, when the central mass reaches 0.38 Mq or 38% 
of the initial core mass). In the aligned case (Model A), the inner accretion flow remains 
dominated by the highly dynamic DEMS at late times, with no sign of RSD formation. 

For the intermediate tilt angle case of 9q = 45° (Model D), the inner structure of the 
protosteUar accretion flow is shaped by the tussle between RSD and DEMS. Fig. [7] shows 
that, at the plotted time. Model D has several spiral arms that appear to merge into a 
rotating disk. There are, however, at least three low-density "holes" near the center of the 
disk: they are the magnetically dominated DEMS. Movies show that the highly variable 
DEMS are generally conflned close to the center, although they can occasionally expand to 
occupy a large fraction of the disk surface. Overall, the circumstellar structure in Model D 
is more disk-like than DEMS-like. We shall call it a "porous disk," to distinguish it from the 
more flUed-in, more robust disk in the orthogonal Model G. Even though most of the porous 
disk has a rotation speed dominating the infall speed, the infall is highly variable, and often 
supersonic. The rotation speed also often deviates greatly from the Keplerian value. Such 
an erratic disk is much more dynamic than the quiescent disks envisioned around relatively 
mature (e.g.. Class II) YSOs. The intermediate tilt angle case drives a powerful bipolar 
outflow, unlike the orthogonal case, but similar to the aligned case. This is consistent with 
the rate of ang ular momentum removal in creasing!: with decreasing tilt angle (i.e., from 90° 



to 45°; see also lCiardi fc Hennebellell2010l ) 



As the strength of the initial magnetic fleld in the core increases, the DEMS becomes 
more dominant. This is illustrated in the middle column of Fig. [71 where the three cases 
with an intermediate fleld strength corresponding to A = 4.86 (and Agflp = 6.85) are plotted. 
In Model H, where the magnetic and rotation axes are orthogonal, a relatively small (with 
radius of ^ 10^ AU) rotationally dominated disk is clearly present at the time shown. As in 
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the weaker field case of Model G, it is fed by prominent "pseudo-spirals" which are part of 
a magnetically-induced curtain in 3D (see the second panel of Fig. [2j). Compared to Model 
G, the curtain here is curved to a lesser degree, which is not surprising because the rotation 
is slower due to a more efficient braking and the stronger magnetic field embedded in the 
curtain is harder to bend. The disk is also smaller, less dense, and more dynamic. It is more 
affected by DEMS, which occasionally disrupt the disk, although it always reforms after 
disruption. Overall, the circumstellar structure in Model H is more RSD-like than DEMS- 
like. As in Model D, we classify it as a "porous disk." As the tilt angle decreases from 
90° to 45° (Model E) and further to 0° (Model B), the rotationally dominated circumstellar 
structure largely disappears; it is replaced by DEMS-dominated structures. Even though 
there is still a significant amount of rotation in the accretion fiow, a dense coherent disk is 
absent. We conclude that for a moderately strong magnetic field of A = 4.86 the formation 
of RSD is suppressed if the tilt angle is moderate. 

In the cases of the strongest magnetic field corresponding to A = 2.92 (and Aeff = 4.12), 
the formation of RSD is suppressed regardless of the tilt angle, as can be seen from the 
last column of Fig. [71 For the orthogonal case (Model I), the prominent "pseudo-spirals" in 
the weaker field cases of Model G (A = 9.72) and H (A = 4.86) are replaced by two arms 
that are only slightly bent. They are part of a well-defined pseudodisk that happens to lie 
roughly in the ^ = (or x-z) plane (see the second panel of Fig. [8]). In the absence of any 
initial rotation, one would expect the pseudodisk to form perpendicular to the initial field 
direction along the x-axis, i.e., in the x = (or y-z) plane. Over the entire course of core 
evolution and collapse, the rotation has rotated the expected plane of the pseudodisk by 
nearly 90°. Nevertheless, at the time shown (when the central mass reaches 0.2 Mq), there is 
apparently little rotation left inside 10^ AU to warp the pseudodisk significantly. Except for 
the orientation, this pseudodisk looks remarkably similar to the familiar one in the aligned 
case (Model C; see the first panel of Fig. [8]). In particular, there are low-density "holes" 
in the inner part of both pseudodisks which are threaded by intense magnetic fields and 
surrounded by dense filaments: they are the DEMS. In the intermediate tilt angle case of 
45° (Model F, not shown in the 3D figure), the pseudodisk is somewhat more warped than 
the two other cases, and its inner part is again dominated by DEMS. It is clear that for a 
magnetic field of A of a few, the inner circumstellar structure is dominated by the magnetic 
field, with rotation playing a relatively minor role; the RSD remains suppressed despite the 
misalignment. 

To better estimate the boundary between the cores that produce RSDs and those that 
do not, we carried out two additional simulations with Oq = 90° (Models P and Q in Table 1). 
We found a porous disk in Model P (A = 4.03 and Agfr = 5.68), as in the weaker field case of 
Model H, but no disk in Model Q (A = 3.44 and Agflp = 4.85), as in the stronger field case of 
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Model 1. From this, we infer that the boundary hes approximately at A 3.5 (or Agff ^ 5). 



Discussion 



5.1. Comparison with Ijoos et al 



Our most important quahtative result is that the misalignment between the magnetic 
field and rotation axis tends to promote the formation of rotationally supported disks, espe- 
cially in we akly magne t ized d ense cores. This is in agreement with the conclusion previously 
reached bv Ijoos et al.l (120121 : JHC12 hereafter), using a different numerical code and some- 
what different problem setup. Their calculations were carried out using an adaptive mesh 
refinement (AMR) code in the Cartesian coordinate system, with the central object treated 
using a stiffening of the equation of state, whereas ours were done using a fixed mesh re- 
finement (FMR) code in the spherical coordinate system, with an effective sink particle 
at the origin. Despite the differences, these two distinct sets of calculations yield qualita- 
tively similar results. The case for the misalignment promoting disk formation is therefore 
strengthened. 

Quantitatively, there appears to be a significant discrepancy between our results and 
theirs. According to their Fig. 14, a Keplerian disk is formed in the relatively strongly 
magnetized case of A = 3 if the misalignment angle Oq = 90°. Formally, this case corresponds 
roughly to our Model I (A = 2.92 and Oq = 90°), for which we can rule out the formation of 
a rotationally supported disk with confidence (see Fig. [Tj). 

We believe that the discrepancy comes mostly from the initial density profile adopted, 
which a ffects the de gree of magnetization near the core center for a given global mass-to-fiux 
ratio A. Ijoos et al.l adopted a centrally condensed initial mass distribution 

Po 



p(r) 



1 + (r/ro)^ 



(8) 



with the characteristic radius r p set to 1 /3 of the core rad ius i?c, so that the central-to-edge 
density contrast is 10 (see also lCiardi fc Hennebelldl2010l ). It is easy to show that, for this 
density profile and a uniform magnetic field, the mass-to-fiux ratio for the fiux tube passing 
through the origin is 



Ae = 



Xq arctana;o 



-A 



(9) 



2(xo — arctanxo) 

where Xq = Rc/tq and A is the global mass-to-fiux ratio for the core as a whole. For xq = 3, 
we have Ac = 3.21 A. In other words, the central part of their core is substantially less 
magnetized relative to mass than the core as a whole, due to the initial mass condensation. 
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For the A = 3 case under consideration, we have Ac = 9.63, which makes the material on 
the central flux tube rather weakly magnetized (relative to mass) . The magnetization of the 
central region is important, because the central par t is accreted first and the star formation 



efl&ciency in a core may not be 100% eflicient (e.g., lAlves et al.l 120071 ). If we deflne as in ^ 
an effective mass-to-flux ratio for the (cylindrical) magnetic flux surface that encloses 1/3 of 
the core mass, then Agff = 2.39A for the density distribution adopted by JHC12 (equation 
jH]). It is signiflcantly different from the effective mass-to-flux ratio of Agff = 1.41 A for the 
uniform density that we adopted. 

The above difference in Agff makes our strongest fleld case of A = 2.92 (Agff = 4.12) more 
directly comparable to JHC12's A = 2 (Agff = 4.78) case. There is agreement that, in both 
cases, the formation of a RSD is suppressed, even when the tilt angle 9o = 90°. These results 
suggest that RSD formation is suppressed when the effective mass-to-flux ratio Agff ^ 5, 
independent of the degree of fleld-rotation misalignment, consistent with the conclusion we 
reached toward the end of ^ based on Models P and Q. 

Similarly, JHC12's A = 3 (Agff = 7.17) models may be more directly comparable to our 
A = 4.86 (Aeff = 6.85) models. Our calculations show that a (more or less) rotationally 
supported disk is formed in the extreme 9q = 90° case for Agff = 6.85 (Model H, see Fig. 
[7j) but not in the intermediate tilt angle 9q = 45° case (Model E). This is consistent with 
their Fig. 14, where a Keplerian disk is formed for = 90° for Agff = 7.17, but not for 
^0 = 45°. In the latter case, JHC12 found a disk-like structure with a flat rather than 
Keplerian rotation proflle. They attributed the flat rotation to additional support from the 
magnetic energy, which dominates the kinetic energy at small radii. We believe that their 
high magnetic energy comes from strongly magnetized, low-density lobes (i.e., DEMS), which 
are clearly visible in our Model E (see Fig. (T]). Although there is still signiflcant rotation on 
the 10^ AU scale in Model E, we flnd the morphology and kinematics of the circumstellar 
structure too disorganized to be called a "disk." We conclude that for a moderate fleld 
strength corresponding to Agff ^ 6-7, a rotationally supported disk does not form except 
when the magnetic fleld is tilted nearly orthogonal to the rotation axis. 

Our weak-fleld case of A = 9.72 (Aeff = 13.7) can be compared with JHC12's A = 5 
(Aeff = 12.0) case. In both cases, a well-deflned rotationally supported disk is formed when 
9o = 90° (see our Fig. Hand their Fig. 12). For the intermediate tilt angle 9q = 45°, JHC12 
obtained a disk-like structure with a flat rotation curve, with the magnetic energy dominating 
the kinetic energy at small radii. This is broadly consistent with our intermediate tilt angle 
Model D, where a highly variable, "porous" disk is formed, with the central part often 
dominated by strongly magnetized, low-density lobes. There is also agreement that the disk 
formation is suppressed if = even for such a weakly magnetized case. It therefore appears 
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that, for the weak-field case of Agff ^ 10, a rotationally supported disk can be induced by a 
relatively moderate tilt angle of 9o > 45°. 

The result that a misalignment between the magnetic field an d rotation axis helps disk 



formation by weakening magnetic braking may be counter-intuitive. iMouschovias fc Paleologou 



( 119791 ) showed analytically that, for a uniform rotating cylinder embedded in a uniform static 
external medium, magnetic braking is much more efficient in the orthogonal case than in 
the aligned case. This analytic result may not be directly applicable to a collapsing core, 
however. As emphasized by JHC12, the collapse drags the initially uniform, rotation-aligned 
magnetic field into a configuration that fans out radially. JHC12 estimated analytically that 
the collapse-induced field fan-out could in principle make the magnetic braking in the aligned 
case more efficient than in the orthogonal case. The analytical estimate did not, however, 
take into account of the angular momentum removal by outffow, which, as we have shown in 
§3.2[ is a key difference between the weak-field (A = 9.72, Agff = 13.7) aligned case (Model 



A) where disk formation is suppressed and its orthogonal count erpart (Model G) that doe s 



produce a rotationally supported disk (see Figs. [T] and [U and also lCiardi fc Hennebelldl2010l ). 



The generation of a powerful outfiow in the weak-field aligned case (Model A) is facil- 
itated by the orientation of its pseudodisk, which is perpendicular to the rotation axis (see 
Fig. [2j). This configuration is conducive to both the pseudodisk winding up the field lines 
and the wound-up field escaping above and below the pseudodisk, which drives a bipolar 
outflow. When the magnetic field is tilted by 90° away from the rotation axis, the pseudodisk 
is warped by rotation into a snail-shaped curtain that is unfavorable to outfiow driving (see 
Fig. [2j). The outfiow makes it more diflicult to form a rotationally supported disk in the 
aligned case. Disk formation is further hindered by magnetically dominated, low-density 
lobes (DEMS), which affect the inner part of the accretion ffow of the aligned case more 
than that of the perpendicular case, at least when the ffeld is relatively weak (see Fig. [1] 
and[7j). For more strongly magnetized cases, the DEMS becomes more dynamically impor- 
tant close to the central object, independent of the tilt angle Oq (see Fig. [8j). DEMS-like 
structures were also seen in some runs of JHC12 (e.g., the case of A = 2 and 9o = 0°; see 



thei r Fig. 19) but were not com mented upon. As stressed previously by IZhao et al.l ( 120111 ) 



and iKrasnopolsky et al.l ( 120121 ) and conffrmed by our calculations, the DEMS presents a 



formidable obstacle to the formation and survival of a rotationally supported disk. 



5.2. Misalignment not Enough for Disk Formation in General 



While there is agreement between JHC12 and our calculations that misalignment be- 
tween the magnetic ffeld and rotation axis is beneffcial to disk formation, it is unlikely that 
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the misalignment alone can enab le disk formation around th e majority of young stellar ob- 
jects. The reason is the following. iTroland fc Crutcherl ( 120081 ) obtained a mean dimensionless 
mass-to-flux ratio of Aeff ^ 2 through OH Zeeman observations for a sample of dense cores 
in nearby dark cloudq^. For such a small Agff , disk f ormation is completely suppressed, even 
for the case of maximum misalignment of 9q = 90°. ICrutcher et al.l ( 120101 ) argued, however, 
that there is a flat distribution of the total fleld strength in dense cores, from Btot ^ to 
some maximum value S^ax ^ 30 /xG; the latter corresponds to Agfr 1, so that the mean Aeff 
stays around 2. If this is the case, some cores could be much more weakly magnetized than 
others, and disks could form preferentially in these cores. However, to form a rotationally 
supported disk, the core material must have (1) an eflFective mass-to-flux ratio Agff greater 
than about 5 and (2) a rather large tilt angle (see discussion in the preceding section). If one 
assumes that the core-to-core variati on of Apff comes mostly from the field strength rather 
than the column density (as done in ICrutcher et al.ll2010l ). then the probability of a core 
having Aeff ^ 5 (or Btot < 6/xG) is ^ 1/5. Since a large tilt angle of Oq ^ 90° is required to 
form a RSD for Aeff ^ 6-7, the chance of disk formation is reduced from ^ 1/5 by at least 
a factor of 2 (assuming a random orientation of the magnetic field relative to the rotation 
axis), to ^ 10% or less. 

The above estimate is necessarily rough, and can easily be off by a factor of two in either 
direction. It is, however, highly unlikely for the majority of the cores to simultaneously satisfy 
the conditions on both Aeff and for disk formation. The condition Aeff ^ 5 is especially 
diflicult to satisfy because, as noted earlier, it implies that the dense cores probed by OH 
observations must have a total fleld strength ^tot 5, 6/iG, less t han o r comparable to the 
well-deflned median fleld strength inferred by iHeiles fc Trolandl (120051 ) for the much more 
diffuse, cold neutral atomic medium (CNM). It is hard to imagine a reasonable scenario in 
which the majority of dense cores have magnetic flelds weaker than the CNM. 



We note that iKrumholz et al.l (120131 ) independently estimated a range of ^ 10-50% 
for the fraction of dense cores that would produce a Keplerian disk based on Fig. 14 of 
JHC12. Their lower limit of ^ 10% is in agreement with our estimate; in both cases, 
the fraction is dominated by weakly magnetized cores that have mass-to-flux ratios greater 
than 5 and moderately large tilt angles. Their upper limit of ^ 50% is much higher than 
our estimate, mainly because it includes rather strongly magnetized cores with mass-to-flux 
ratios as small as 2. Since our calculations show that such strongly magnetized cores do not 



^We assume that the mass-to-flux ratio measured by Troland & Crutcher on the core scale is the same 
as the effective mass-to-flux ratio near the core center because, for significantly magnetically supercritical 
cores, ambipolar diffusion is generally ineffective in reducing the value of A near the center relative to that 
in the envelope. 
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produce rotationally supported disks even for large tilt angles, we believe that this upper 
limit may be overly generous. 

Whether dense cores have large tilt angle bet ween the magnet ic field and rotation 
axis that are conducive to disk formation is unclear. IHuU et al.l ( 120121 ) measured the field 
orientation on the 10^ AU-scale for a sample of 16 sources using millimeter interferometer 
CARMA. They found that the field orientation is not tightly correlated with the outfiow 
axis; indeed, the angle between the two is consistent with being random. If the outfiow 
axis is aligned with the core rotation axis and if the field orientation is the same on the 
core scale as on the smaller, 10^ AU-scale, then Oq would be randomly distributed between 
and 90°, with half of the sources having 9q > 60°. However, the outfiow axis may not 
be representative of the core rotation axis. This is because the (fast) outfiow is thought to 
be driven m aRnetocentrifuRally from the inner part of the circumstellar disk (on the AU- 
scale or less; IShu et al.ll2000l : iKonigl fc Pudritzll2000l ). A parcel of core material would have 
lost most of its angular momentum on the way to the outfiow launching location; the torque 
(most likely magnetic or gravitational) that removes the angular momentum may also change 
the direction of the rotation axis. Similarly, the field orientation on the 10^ AU-scale may 
not be representative of the initial field orientation on the larger core scale. The magnetic 
field on the 10^ AU scale i s more prone to distor tion by collapse and rotation than that 
on the core scale. Indeed, IChapman et al.l ( 120131 ) found that the field orientation on the 
core scale measured using single-disk telescope CSO is within 20° of the outfiow axis for 3 
of the 4 sources in their sample; the larger angle measured in the remain source may be 
due projection e fiFects because its outfiow axis is close to the line-of-sight. If the result of 
Chapman et al.l is valid in general and if the outfiow axis refiects the core rotation axis, then 
dense cores with large tilt angle would be rare. In this case, disk formation would be 
rare according to the calculations presented in this paper and in JHC12, even in the unlikely 
event that the majority of dense cores are as weakly magnetized as Agff ^ 5. 

Rotat ionally supported disks are observed, however, routinely around evolved Class II 
YSOs ("seelWilliams fc Ciezall201llfor a recent review), and increasingly around y ounRer Class 



I (e.g . Jj^rgensen et al. 



2009 



Le 



201 ll : Takakuwa et al. 2012) and even Class ( ITobin et al. 



20121) sources. When and how such disks form remain unclear. Our calculations indicated 
that the formation of large observable (10^ AU-scale) rotationally supported disk^ is difiicult 
even in the presence of a large tilt angle during the early protostellar accretion (Class 0) 



^We cannot rule out the existence of small, AU-scale disks, because the size of our effective "sink particle" 
is 6.7 AU. Small disks n iay be needed to drive fast outflows during the Class phas e, and may form through 
non-ideal MHD effects (|Machida et al.lboioi [Papp Basulboioi IPapp et al.ll2Q12l ). 
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phasqj. This may not contradict the available observations yet ( IMaury et al.ll2010). since 



only one Keplerian disk is found around a (late) Class source so far (L1527; 



Tobin et al. 



20121); it could result from a rare combination of an unusually weak magnetic field and a 



large tilt angle ^o- If it turns out, through future observations using ALMA and JVLA, 
that large-scale Keplerian disks are prevalent around Class sources, then some crucial 
ingredients must be missing from the current calculations. Possible candidates include non- 
ideal MHD eflFects and turbulence. The existing calculations indicate that realistic levels of 



large-scale disk formation {] 


VIellon & Li 


2009; 


Li et al. 


2011; 


KrasnoDolskv et al. 


2012 


; see 


also 


Krasnopolskv & Konipl 


2002|and 


Braidinp & Wardle 


2012a, 


q), although misalignment 



has yet to be considered i n such calculations. Sup e rsonic turbulence was found to be con- 



ducive to disk format ion (ISantos-Lima et al.l 120121 . l2013l : ISeifried et al.l l2012l : iMyers et al. 



20121; IJoos et al.ll2013l ). although dense cores of low-mass star formation typically have sub- 



sonic non-thermal line-width and it is unclear whether subsonic turbulence can enable disk 
formation in dense cores magnetized to a realistic level. If, on the other hand, it turns 
out that large-scale Keplerian disks are rare among Class sources, then the question of 
disk growth becomes paramount: how do the mostly undetectable Class disks become 
detectable in the Class I and II phase? If the magnetic braking plays a role in keeping the 
early disk undetectable, then its weakening at later times may promote rapid disk growth. 
One possibility for the late weakening of magnetic brakin g is the deplet ion of the protostella r 
envelope, either by outfiow stripping ( IMellon fc Lill2008l ) or accretion (IMachida et al.ll201ll ). 
It deserves to be better quantified. 



5.3. Summary 



We carried out a set of MHD simulations of star formation in dense cores magnetized to 
different degrees and with different tilt angle s between the ma gnetic field and the rotation 
axis. We confirmed the qualitative result of Ijoos et al.l ( 120121 ) that misalignment between 
the magnetic field and rotation axis tends to weaken magnetic braking and is thus conducive 
to disk formation. Quantitatively, we found however that the misalignment enables the 
formation of a rotationally supported disk only in dense cores where the star- forming material 
is rather weakly magnetized, with a dimensionless mass-to-fiux ratio > 5; large misalignment 
in such cores allows the rotation to wrap the equatorial pseudodisk in the aligned case into 



^We are unable to follow the collapse until most of the core material is accreted into the central object for 
numerical reasons. Such a complete acc retion may not be realistic, however, if the star formation efficiency 
in a core is typically as low as 1/3 (e.g. jAlves et al.ll2QQ7[ ). 
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a curved curtain that hinders outflow driving and angular momentum removal, making 
disk formation easier. In more strongly magnetized cores, disk formation is suppressed 
independent of the misalignment angle, because the inner part of the protostellar accretion 
flow is dominated by stron gly magnetized, low-density regions. If dense cores are as strongly 



magnetized as inferred by iTroland &: Crutcherl (120081 : with a mean mass-to-flux ratio ~ 2) 



it would be diflicult for the misalignment alone to enable disk formation in the majority of 
them. We conclude that how protostellar disks form remains an open question. 

We thank Mark Krumholz for useful discussion. This work is supported in part by NASA 
grant NNX10AH30G and the Theoretical Institute for Advanced Research in Astrophysics 
(TIARA) in Taiwan through the CHARMS project. 
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Fig. 5. — Infall (left panel) and rotation (middle) speed and the specific angular momentum 
(right) on the equatorial plane along 4 representative directions for the aligned (dashed 
lines) and orthogonal (solid) case. For comparison, the free-fall speed for the aligned case 
and the Keplerian profile for the orthogonal case are plotted as dotted-dash line in the left 
and middle panel respectively. It is evident that a rotationally supported disk is formed in 
the orthogonal but not the aligned case. In the latter case, the specific angular momentum 
of the equatorial flow decreases significantly upon entering the pseudodisk at ^ 10^^ cm and 
the DEMS-dominated region around ^ 10^^ cm. 



Fig. 6. — Left: Same as Fig. [3l except that the color map is for the toroidal component 
of the equatorial magnetic field (rather than density), showing the inward spiraling of two 
strains of magnetic field lines of opposite polarity. For clarify, the plotted field strength is 
capped from above and below at ±3 x 10~^G. Right: The plasma /3 along 4 representative 
directions in the equatorial plane showing that the inner 100 AU region in the orthogonal 
case (i.e., on the Keplerian disk; solid lines) is dominated by thermal rather than magnetic 
pressure, whereas the opposite is true for the aligned case (dashed). 
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Fig. 7. — Snapshots of Models A~I (see Table 1 for model parameters) when the central mass 
reaches 0.2 Mq. Plotted are maps of the logarithm of density (color range: —18 to —12) and 
velocity unit vectors on the equatorial plane. The box size in each panel is 10^^ cm. 
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Fig. 8. — Same as in Fig. O but for the stronger field (A = 2.92) cases of Model C (aligned) 
and Model I (orthogonal) at a time when the central mass reaches 0.2 M©. For clarity, 
only field lines originated from the bottom x-y and left y-z plane are plotted, respectively. 
Note the difference in orientation of the pseudodisks in the two cases. The inner part of 
both pseudodisks is dominated by low-density, strongly magnetized regions. The length is 
in units of AU. 



